Magnetoelastic coupling in BaTiOa-based multiferroic structures 



OO 

o 
o 

> 
O 



o 



> 



00 

o 



X 



C. A. F. V&E 

Department of Applied Physics, Yale University, New Haven, Connecticut 06520 and 
Center for Research on Interface Structures and Phenomena (CRISP), 
Yale University, New Haven, Connecticut 06520 
(Dated: November 13, 2008) 

Analytical expressions for the magnetoelastic anisotropy constants of cubic magnetic systems 
are derived for rectangular and oblique distortions originating from epitaxial growth on substrates 
with lower crystal symmetry. In particular, the temperature variation of the magnetic properties 
of magnetic films grown on barium titanate (BaTiOs) substrates are explained in terms of strain- 
induced magnetic anisotropics caused by the temperature dependent phase transitions of BaTiOs. 
Our results quantify the experimental observations in ferromagnet/BaTiOs-based structures, which 
have been proposed as templates for magnetoelectric composite heterostructures. 
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I. INTRODUCTION 

Complex oxides are characterized by a strong coupling 
between charge, spin and lattice distortions, giving rise to 
a rich variety of electronic and magnetic behavior which 
is at the origin of the multifunctional properties exhib- 
ited by this class of materials. In particular, much effort 
is being made towards the control of these properties by 
means of electric and magnetic fields. While intrinsic 
multiferroic materials already exhibit a coupling between 
magnetic and electric order parameters,-^'^'^ the magni- 
tude of the coupling in these compounds is weak, and a 
new class of composite materials is being developed that 
exploit the different properties of the constituent com- 
pounds to tailor systems with larger susceptibilities and 
more varied response functions than those of single com- 
ponent systems. 

In this context, the cross coupling in ferroelec- 
tric/magnetic heterostructures and composite materials 
is being intensively studied. The coupling mechanism 
can be classified as strain-driven or charge driven: in 
the first case, the piezoelectric effect of ferroelectrics is 
used to modify the magnetic properties of ferro- or fer- 
rimagnets via magnetoelastic coupling;'''^'^'^'^ in the sec- 
ond case, the coupling originates from modifications in 
the electronic properties of the magnetic system driven 
by charge modulation induced by the polarization state 
of the ferroelectric.^ 

A particular case of strain-mediated magnetoelectric 
coupling occurs in epitaxial heterostructures composed 
of thin magnetic films grown on BaTiOa single crystals, 
including l,&^_^SvM^Ozr^^^^^ SrRuOS^ CoFe204?i^ 
Fej^ and Ye^O/i^i^^^ where the changes in strain oc- 
cur as a function of temperature via the phase transi- 
tions of BaTiOa or via applied electric fields. However, 
to date, the magnetoelastic coupling has not been ad- 
dressed explicitly in detail. Here, we present analytical 
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expressions for the magnetoelastic contributions to the 
magnetic anisotropy induced by elastic coupling between 
the different phases of the BaTiOs substrate to the epi- 
taxial magnetic film. This coupling explains the mag- 
netic behavior exhibited by these systems as a function 
of temperature in terms of the changes in the magnetoe- 
lastic anisotropics. 

BaTiOs has a simple perovskite structure at high tem- 
peratures and undergoes several crystal phase transitions 
as a function of temperatureiiS'i^ It is cubic in the tem- 
perature range from ^1733 K to ^393 K (a = 3.996 
at 393 K), tetragonal at temperatures down to ~278 K 
(a = 3.9920 A and c = 4.0361 A at 293 K), orthorhom- 
bic down to ~183 K {a^ = 5.669 A, bo = 5.682 A, and 
Co = 3.990 A at T = 263 K), and rhombohedral at 
lower temperatures (a = 4.001 A, 7 = 89.85° at 105 K). 
For temperatures below ~393 K, BaTiOs is ferroelectric. 
The non-cubic lattices correspond to slight distortions 
of the cubic phase, with the direction of the distortion 
linked to that of the electric polarization, p-,i2iiiS in the 
tetragonal phase the electric polarization points along 
the c axis; in the orthorhombic phase, the spontaneous 
polarization points along the (110) directions of the cu- 
bic cell, leading to a slight elongation of the unit cell 
along that direction. Although the symmetry group is or- 
thorhombic, with lattice parameters Oo, bo, and Co when 
P points along [110], it is often convenient to consider 
the monoclinic cell with = = (a^ -|- b'^Y^'^ /2 and 
7 = 2 arctan(ao/6o), giving a = b = 4.013 A, c = 3.99 A, 
7 = 89.869° at T = 263 K. For the rhombohedral phase, 
the polarization points along the (111) directions of the 
cubic cell, corresponding to an elongation of the cubic 
cell along P. 



II. MAGNETO-ELASTIC ENERGY 
EXPRESSIONS 

In the following, we are interested in determining the 
effect of the in-plane strain on the anisotropy constants 
of ferromagnetic (or ferrimagnetic) films due to magne- 
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FIG. 1: Schematic of a general (oblique) distortion of a rect- 
angular unit cell. The strain components correspond to the 
normalized difference between the distorted and the original 
equilibrium unit cell. 



toelastic coupling. We assume that the strain propa- 
gates uniformly throughout the film. When relaxation 
occurs, as is often the case in real systems, the expres- 
sions below correspond to upper bounds to the magneto- 
lastic anisotropy induced in the magnetic system. If the 
residual strain is known, the effective anisotropics can be 
estimated. 

To calculate the magnetoelastic anisotropy contribu- 
tions arising from strain, we minimize the sum of the 
elastic (Fdas) and magnetoelastic (Fi„o) energies for a 
given strain state of the system. We consider the case of 
a cubic magnetic system>2£i2iiS^ 
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where at are the direction cosines of the magnetization, 
eij are the components of the strain tensor, bi are the 
magnetoelastic coupling coefficients, Cy are the elastic 
constants and i,j are the cartesian coordinate indexes. 
The magnetoelastic constants are related to the magne- 
tostriction coefficients byi^ 



h = -3(cii - ci2)Aioo/2 
62 = — 3C44A111 



(3) 
(4) 



For a general in-planc distortion in a cubic lattice, it 
is straightforward to calculate the in-plane strain tensor 
components, since by definition the latter correspond to 
the relative deviations from the equilibrium lattice con- 
stant; with the aid of Fig. [2 we find: 

e^x = (acosCi ~ ao)/ao (5) 
eyy = (CC0SC2 - co)/co (6) 
exy = [(a/ao)sinCi + (c/co)sinC2]/2, (7) 

where qq, Cq are the equilibrium in-plane lattice param- 
eters of the magnetic film and a, c, 7 and are the 



parameters characterizing the planar lattice distortion. 
Here, we discuss the case of a (001) ferromagnetic film, 
where — cq, but the calculations can be easily modi- 
fied for other cases. The out of plane strain components 
are free to settle to the values that minimize the total 
energy (zero out of plane stress). The magnetoelastic 
anisotropy is then calculated by finding the minimum of 
-Fm = -Fme + -Feias with respect to exz, eyz and e^z- 
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inserting ((8] fT0|) into the expression for i^m gives, retain- 
ing only the terms that depend on the magnetization: 
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where the last term is strain independent, correspond- 
ing to a sum of terms of the order of bl/cu and 62/^44 
with a1 and dependencies. These are much smaller 
than the first terms and will be ignored in the follow- 
ing discussion. They arise from a shear strain exerted 
by the magnetization on the lattice (which favors spe- 
cific directions that depend on the relative sign of the 
magnetoelastic coupling coefficients) and are usually ne- 
glected to first order of approximation. We note that the 
magnetoelastic terms are uniaxial and lower in symmetry 
than the magnetocrystalline anisotropy of the unstrained 
cubic magnetic system. 

A more convenient expression is obtained by using 
spherical coordinates, where the polar axis is set along 
the out of plane direction: 



[bi{exx - £yy) cos 2(j) + b2£xy sm 2(j)\ sin^ 6 - 



bi 



(l + 2ci2/cn)(e^ 



ej^j^)cos 9 
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where the first two terms correspond to in-plane magne- 
toelastic anisotropics, and the last term is the perpendic- 
ular magnetoelastic anisotropy contribution. An in-plane 
uniaxial anisotropy arises for uniaxial in-plane strains 
i^xx 7^ £yy) or for shear strains, although the anisotropy 
direction is distinct for these two cases: in the first case, 
it is oriented along the main crystal axes, while for shear 
strains, the anisotropy is directed along the in-plane diag- 
onal directions. A perpendicular anisotropy arises when 
fizz 7^ 0. The relative magnitude of biezz with respect to 
the magnetostatic energy term (shape anisotropy, which 
favors in-plane magnetization), which is on the order of 
27rM^, where is the saturation magnetization, de- 
termines whether a state of in-plane or perpendicular 
magnetization is favored. For the particular case of a 
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biaxial strain (e^x = Cyy, £xy = 0), the magnetoelastic 
anisotropy energy reduces to a perpendicular magnetic 
anisotropy term only, with Kp — — foi(l + 2ci2/cii)exx, 
as expected.^ 

We consider next the effect of square, rectangu- 
lar, oblique and rhomboidal distortions on the mag- 
netoelastic anisotropy. These planar domains arise 
in the cubic, tetragonal, orthorhombic and rhombo- 
hedral phases of BaTiOs, and are responsible for 
the large magnetic anisotropy modulations observed in 
ferromagnetic/BaTiO3(001) heterostructures. We con- 
sider the magnetoelastic anisotropy constants, Kt, Kq, 
Kp, expressed as: 

= [Kt cos 20 + Ko sin 2(/)] sin^ + Kp cos^ 0. (13) 

The direction of easy magnetization depends on the par- 
ticular K values; positive Kp favors in-plane magnetiza- 
tion. The magnetoelastic anisotropy for the cubic phase 
of BaTiOs {(i = C2 = 0, a = c) leads to zero in-plane 
anisotropy. For the tetragonal phase (Ci = C2 — 0), 
two types of domains are possible, c-domains (square) 
and a-domains (rectangular, which may co-exist in two 
orthogonal directions). A multidomain state can arise 
when cooling the BaTiOs from the cubic to the tetrag- 
onal phase after film growth at elevated temperatures, 
for example. For this phase, no shear strains are present 
and Ko = 0. The orthorhombic phase may give rise to 
(two) rectangular (Ci = ^2 = 0, a ^ c) and (two) oblique 
(Ci — (2, a — c) domains. Finally, the rhombohedral 
phase (Ci = C21 a = c) is characterized by (two) oblique 
domains. It also follows that for the oblique domains of 
the orthorhombic and rhombohedral phases of BaTiOa, 
Kt^O. 



III. APPLICATION TO EPITAXIAL 
MAGNETIC FILMS 

In general, it may be difficult to control the ferroelec- 
tric state of BaTiOa across the different phase transi- 
tions. For example, the growth of oxide magnetic films 
is often carried out at elevated temperatures, where 
BaTiOa is cubic and, as the temperature is lowered to 
room temperature, the BaTiOa tends to break into a 
ferroelectric multidomain state, giving rise to a and c 
domains in a BaTiO3(001) single crystal. Similar pro- 
cesses operate at the other phase transitions, leading to 
the presence of different types of planar domains at the 
interface with the magnetic film. In such cases, an effec- 
tive magnetoelastic anisotropy corresponds to an average 
over the relative number and size of domains present; 
in the case where an equal number of the two possible 
oblique domains are present in the rhombohedral phase, 
no effective in-plane anisotropy would result. However, 
experimental estimates of the relative population of the 
different domains is difficult. In the following we will 
consider single domain states; for similar, but orthogo- 
nal, domains, the in-plane anisotropy constants will be 



identical, but with opposite signs, while the perpendic- 
ular anisotropies are unaffected. Another difficulty re- 
lates to the fact that for systems with large lattice mis- 
match, such as Fe304 and CoFe204 grown on BaTiOa, 
the strain is strongly reduced for films thicker than the 
coherent critical thickness through the onset of misfit 
dislocations;-2M5 in the case of Fe3O4/BaTiO3(100), the 
experimental evidence points to full relaxation at the 
temperature of growthJi Here, we do not take explicitly 
into account the effect of dislocations to the magnetoe- 
lastic energy, but instead account for strain relaxation by 
adjusting the lattice parameter of the magnetic film so 
that the correct strain value is obtained at the measure- 
ment temperature, as will be discussed in more detail 
below. We consider three magnetic systems that have 
been studied experimentally and which are of particular 
interest: CoFe204,i^ FesOi^'^ and Fe films^'' grown 
on BaTi03. In Table HI we list the relevant materials 
parameters for these magnetic systems. 



A. Fe3O4/BaTiO3(100) epitaxial films 

We consider first the case of Fe3O4/BaTiO3(100) epi- 
taxial films, where experimentally it is found that the 
Fe304 has a residual strain at room temperature with 
a sign that is opposite that expected from the lattice 
mismatch with BaTi03; this has been explained as due 
to strain relaxation at the temperature of growth, fol- 
lowed by an in-plane dilation due to the abrupt increase 
in the cell parameters of BaTiOs at the cubic to tetrag- 
onal phase transition.^- Therefore, the magnetoelastic 
anisotropy contribution to the magnetic energy is very 
different from that expected for a fully strained Fe304 
film: while the values for the magnetostriction coef- 
ficients of Fe304 favor a state of in-plane magnetiza- 
tion, experimentally a non-zero perpendicular remanence 
and magnetic anisotropy are observedJ^iii In Fig. [2] we 
present our estimate for the variation of the magnetoe- 
lastic anisotropy constants as a function of temperature, 
where we took into account the temperature variation 
of the lattice constants of BaTi03 and FesO^Mi^ The 
strain relaxation was taken into account, in a first ap- 
proximation, by shifting the lattice constant of magnetite 
such that the out of plane strain at room temperature co- 
incides with the experimental value of -0.62% for tetrag- 
onal domains.^'' This assumes that the misfit dislocation 
density remains constant throughout the temperature 
range and across the phase transitions of BaTiOs. One 
indication that this approach is a simplification over the 
real processes occurring in the magnetite film is the fact 
that, as observed in Fig. [21 the strain in the cubic phase 
does not drop to zero as the temperature approaches the 
growth temperature, indicating that a strong tempera- 
ture dependence of the lattice constant, or misfit dislo- 
cation density, may occur in the high temperature range. 
The temperature variation of the magnetic anisotropies is 
instructive, however, and helps explain the variations in 
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TABLE I: Room temperature material parameters for Fe, Fe304 and CoFe204. The elastic constants ctj are given in erg/cm'^, 
the magnetization Ms in emu/cm"^, and the cubic magnetocrystalline anisotropy Ki in erg/cm^ (a is the lattice constant). 



Material 


a (A) 






10-^^C44 


lO'^Aioo 


10''Aiii 


Ms 




Fe» 


2.866 


2.322 


1.356 


1.170 


24.1 


-22.7 


1717 


0.48 




8.397 


2.73 


1.06 


0.97 


-19.5 


77.6 


471 


-0.110 


CoFe2 04 


8.381' 


2.57" 


1.500" 


0.853" 


-590'' 


120"^ 


425" 





"From Ref. 26. 
'From Ref. pA 
Trom Ref. |2|. 

"From Ref. for single crystalline Coo.8Fe2.2O4. 
'From Ref. M- 

^From Ref. 23, for sample with composition Coi.iFei,904. 
fFrom Ref. M- 



the magnetic properties observed experimentally. Focus- 
ing on the out of plane anisotropy constant Kp first, we 
see that negative out of plane strains favor perpendicular 
anisotropy, which are large for the rectangular domains 
of the tetragonal phase and for the oblique domains of the 
orthorhombic and rhombohedral phases; they are smaller 
by a factor of about 2 for the square domains of the 
tetragonal phase and the rectangular domains of the or- 
thorhombic phase. Depending on the particular domain 
path followed by the BaTiOa substrate, it is clear that 
large changes in the magnetic anisotropics may follow. 
For example, the experimental results of Vaz et at»^ can 
be explained if, starting from a rectangular domain at 
room temperature, the BaTiOs surface reverts to a rect- 
angular domain in the orthorhombic phase, correspond- 
ing to a drop in the perpendicular magnetic anisotropy, 
followed by an increase in Kp in the rhombohedral phase. 
The in-plane anisotropics are also significant, but they 
tend to be smaller than the perpendicular anisotropics. 
The presence of multidomains in BaTiOs also tends to 
average out the in-plane effective anisotropy. 

It must be emphasized that a negative Kp is not suffi- 
cient to yield perpendicularly magnetized films, since the 
in-plane magnetostatic (shape) anisotropy must also be 
overcome, and the latter is on the order of 2ttM'^ — 1.39 x 
10^ erg/cm^ for magnetite, assuming bulk magnetization 
(thin magnetite films are known to exhibit a slightly re- 
duced magnetization, which would lead to a smaller mag- 



netostatic energy 



33.34.35 



This value would imply that 



the Fe304/BaTi03 films would remain in-plane magne- 
tized, but this assumes a state of uniform strain; in real 
films, one may expect larger interface strains that could 
give rise to larger magnetoelastic anisotropics than those 
predicted in our simple model calculations. 



B. CoFe2O4/BaTiO3(100) epitaxial films 

We consider here the case when the CoFe204 film is 
fully strained to the BaTiOa substrate, a situation which 
is expected to occur at thicknesses below the coherent 
growth thickness. Large changes in magnetization as 
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FIG. 2: Variation of the magnetoelastic anisotropy constants 
of partially strained Fe3O4/BaTiO3(001). Square symbols 
correspond to anisotropy in planar square domains, diamonds 
to oblique domains, and rectangles to rectangular domains; 
full symbols are for in-plane anisotropies, empty symbols for 
the perpendicular anisotropy constants. 



a function of temperature have been reported for par- 
tially relaxed 50-150 nm CoFe2O4/BaTiO3(100) epitaxial 
filmSii^ Cobalt ferrite (CoFe204) is characterized by very 
large magnetostriction coefficients and has been a mate- 
rial of choice for strain-driven composite multiferroicsi^ 
It is ferrimagnetic with a Curie temperature of 790 K and 
is predicted to be nearly fully spin polarizedj ^^'^^ 

The calculated magnetoelastic anisotropies across the 
different phases of BaTiOa using the expressions de- 
rived above are shown in Fig. [3l where we took into ac- 
count the temperature variation of the lattice parameter 
of CoFe204.'^®''^^ The out of plane anisotropies for the 
fully strained film are very large, of the order of 9 x 10^ 
erg/cm'^, and favor in-plane magnetization, so we have 
omitted these from the figure. In the cubic phase, the bi- 
axial strain contributes only to second order to the mag- 
netic anisotropy,'*" and the in-plane magnetic anisotropy 
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will be dominated by the magnetocrystalline anisotropy, 
of the order of 2 x 10^ erg/cm'^ (favoring the in-plane 
(100) directions). In the tetragonal phase, the rectan- 
gular domains induce a very large uniaxial anisotropy 
favoring the direction of in-plane compression [x axis in 
Fig. [1]). In the orthorhombic phase, the distortion in- 
duced by the rectangular domains is smaller, and the 
magnetic anisotropy decreases slightly. The oblique dis- 
tortions of the orthorhombic and rhombohedral phases 
induce relatively small magnetic anisotropics due to the 
much smaller value of Am, and the magnetic anisotropy 
may be dominated by the magnetocrystalline anisotropy, 
assuming that a value comparable to the tabulated room 
temperature value still holds. 




100 200 300 400 500 
Temperature (K) 

FIG. 3: Variation of the in-plane magnetoelastic anisotropy 
constants of fully strained CoFe2 04/BaTi03 (100). Dia- 
monds correspond to anisotropy in oblique domains, rectan- 
gles to rectangular domains; the dashed line corresponds to 
the room temperature magnetocrystalline anisotropy of bulk 
C01.1Fe1.9O4 (which favors the (100) directions). 

Magnetization versus temperature measurements of 
epitaxial CoFe2O4/BaTiO3(001) films have been pre- 
sented by Chopdekar and Suzuki^'^ for a 56 nm film under 
an applied magnetic field iJ = 50 kOe, along the nominal 
[100] and [010] in-plane directions. Even at this very large 
magnetic field, large jumps in the magnetization are ob- 
served at the tetragonal to orthorhombic and orthorhom- 
bic to rhombohedral phase transitions of BaTiOs. As the 
temperature is lowered from 350 K to 150 K, the magne- 
tization for H II [010] direction increases slightly at 290 
K, and then decreases more drastically at 190 K; when 
H II [100] direction, the magnetization behaves in the 
opposite fashion, starting and ending at values that are 
close to those of the magnetization when the magnetic 
field is applied along the [010] direction. These changes 
in the magnetization can be understood by assuming 
predominantly square domains in the tetragonal phase 
of BaTiOa, which give similar anisotropics for both the 



[100] and [010] directions; the presence of a small number 
of rectangular domains, or a small imbalance in the two 
possible rectangular domains, would give rise to a small 
uniaxial magnetic anisotropy favoring, say, the [010] di- 
rection; these domains would either remain rectangular 
in the orthorhombic phase or revert to oblique domains. 
In the orthorhombic phase, the square domains change to 
rectangular and the increase in magnetic energy along the 
[100] direction (leading to a smaller magnetization) could 
result from rectangular domains that are predominantly 
elongated along the [100] direction (the easy axis direc- 
tion is along the direction of compression for CoFe2 04, 
which has bi > 0). The transition to the rhombohedral 
phase would lead to similarly small magnetic anisotropics 
for both [100] and [010] directions, leading to an increase 
in magnetic energy for the [010] direction, and a decrease 
along the [100] direction. Despite clear hysteresis in the 
ferroelectric transitions of BaTiOs, the observation of re- 
producible paths in the magnetic behavior suggests that 
also the ferroelectric domains of BaTiOa also follow sim- 
ilar switching paths with thermal cycling. 



C. Fe/BaTiO3(100) epitaxial films 

As a final example, we consider the case of Fe, where 
the relatively small lattice mismatch with BaTiOa of - 
1.4% should allow the growth of fully strained epitaxial 
films. Further interest in this particular system stems 
from recent theoretical predictions of surface magneto- 
electric effects induced by spin-polarized screening at 
the Fe/BaTiOa interface4i The calculated magnetoelas- 
tic anisotropics for the different phases of BaTiOa are 
shown in Fig. Due to the strong crystal field of metal- 
lic Fe, the orbital moment is quenched, and the magnetic 
anisotropics are relatively small in this system. It can be 
seen that for Fe/BaTiO3(100), the perpendicular mag- 
netoelastic anisotropics favor perpendicular magnetiza- 
tion, but given the very large magnetostatic energy of Fe 
(27rM2 = 1.85 X 10^ erg/cm^), an in-plane magnetized 
state is energetically more favorable. The magnetoelastic 
in-plane anisotropics are also relatively small due to the 
small strains involved, but are likely to strongly influence 
the switching behavior of the Fe film. 

Sahoo et alr^ have reported magnetic measurements 
on a 10 nm polycrystalline Fe film deposited on a 
BaTiO3(100) substrate at 373 K, where the magneti- 
zation under a small applied magnetic field (20-80 Oe) 
shows sudden jumps at the phase transitions of BaTiOa. 
Although the Fe film is likely to be fully relaxed, the 
relative changes in magnetic anisotropy can still be com- 
pared qualitatively with the results shown in Fig. U) For 
the particular field direction, which is applied in the plane 
of the film, the magnetization is found to increase steeply 
with decreasing temperature at both ~278 K and ~183 
K, indicating a decrease in the magnetic anisotropy along 
this direction. The results shown in Fig. 0] suggest that 
the magnetic field is being applied along the a direc- 



6 



0.5 



0.0 

a 

o 

1^-0.5 
o 

^ -1.0 



■1.5 



tetragonal 



romb. 



tDDmnmnnnnnnnnnEDni 

orthor. 



cubic 



an O K 



^□□□□□1 



□□□□□□n°° 



100 200 300 400 
Temperature (K) 



500 



FIG. 4: Variation of the magnetoelastic anisotropy constants 
of fully strained Fe/BaTiO3(100). Square symbols correspond 
to anisotropy in planar square domains, diamonds to oblique 
domains, and rectangles to rectangular domains. Filled sym- 
bols are for in-plane anisotropies and empty symbols are for 
the perpendicular anisotropy constants. 



tion of the rectangular domains of the tetragonal phase, 
which corresponds to a hard axis direction (in Fe this 
is the direction of compression, since bi < 0) and that 
the dominant planar domains in the orthorhombic phase 
of BaTiOa are rectangular, compressed in this direction. 
This also agrees with the measured temperature variation 
of the perpendicular magnetization, where an increase at 
~278 K is followed by a decrease at ~183 K with decreas- 
ing temperature, in agreement with the variation found 



for Kp, shown in Fig. [4] 



IV. CONCLUSIONS 

Explicit expressions for the magnetoelastic anisotropy 
constants have been derived for arbitrary in-plane dis- 
tortions of the unit cell of cubic systems and applied to 
the particular case where those distortions are induced 
by elastic coupling to BaTiOs substrates. BaTiOa un- 
dergoes a series of phase transitions as a function of tem- 
perature and has been employed to study the effect of 
strain on the magnetic properties of epitaxial magnetic 
films and, in particular, aimed at designing new magne- 
toelectric composite systems. Since the elastic proper- 
ties of ferroelectrics, and BaTiOs in particular, can be 
modulated by electric fields via the piezoelectric effect, 
modulation of magnetism via electric fields is possible. 
We have considered three particular magnetic systems 
that have been studied experimentally (Fe304, CoFe204, 
and Fe), and explained the observed magnetic behavior 
in terms of the changes in magnetic anisotropies arising 
from strain induced in the different phases of BaTiOa. 
This work provides a quantitative account of the exper- 
imental observations and opens the way to a more cal- 
ibrated targeting of the properties of heterosystems for 
optimum magnetoelectric coupling effects. 
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